Nonlinear Temporal Dynamics of Strongly Coupled Quantum Dot-Cavity System 
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We theoretically analyze and simulate the temporal dynamics of strongly coupled quantum dot- 
cavity system driven by a resonant laser pulse. We observe the signature of Rabi oscillation in 
the time resolved response of the system (i.e., in the numerically calculated cavity output), derive 
simplified linear and non-linear semi-classical models that approximate well the system's behavior 
in the limits of high and low power drive pulse, and describe the role of quantum coherence in the 
exact dynamics of the system. Finally, we also present experimental data showing the signature of 
the Rabi oscillation in time domain. 
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A single quantum dot (QD) coupled to a photonic crys- 
tal micro-cavity constitutes an integrated nano-photonic 
platform for probing solid state cavity quantum electro- 
dynamic (QED) effects [1]. The eigenstates of this cou- 
pled system form an anharmonic ladder, which results 
in an optical nonlinearity at a single photon level. In 
recent years, this nonlinearity has been used to perform 
all-optical [2, 3] and electro-optic switching [4] as well as 
to generate non-classical states of photons [5-7]. 

In this paper, we study the temporal dynamics of the 
coupled dot-cavity system driven by a short laser pulse 
(Fig. 1 a) using a full quantum optical numerical sim- 
ulation. The oscillatory behavior of the cavity output 
(Fig. 1 b), which is caused by the vacuum Rabi split- 
ting, is analyzed at low, intermediate, and high inten- 
sity of the driving laser. Specifically, we derive a linear 
semi-classical description of the system, and show that 
under weak driving, the coupled QD-cavity system fol- 
lows the same dynamics as a set of two classical linear 
coupled oscillators. Following this, we describe an im- 
proved, non-linear semi-classical model, that mimics the 
quantum optical model very well for both very low and 
high peak intensity of the driving pulse. However, the 
non-linear semi-classical model deviates from quantum 
optical description at intermediate peak intensities of the 
drive pulse and we show that this discrepancy arises from 
the coherence present in the quantum optical system. Fi- 
nally, we present a study of the temporal dynamics as a 
function of the major parameters describing the cavity- 
QD system as well as experimental data showing the sig- 
nature of the Rabi oscillation in time domain. Under 
rotating wave approximation, the quantum-mechanical 
Hamiltonian 1-L describing the coherent dynamics of the 
coupled system is given by 



1-L = cJcjCrV + uOcO^ a + ig{aa^ 



(1) 



Here, Uc and cUa are, respectively, the resonance frequen- 
cies of the cavity and the QD; a is the annihilation oper- 
ator for the cavity mode; a = \g){e\ is the lowering oper- 
ator for the QD with excited state |e) and ground state 



1^); g is the coherent interaction strength between the QD 
and the cavit and h is set to 1. When this system is coher- 
ently driven by a laser pulse with strength Q{t) = ^op{t) 
and a center frequency uji^ the driven Hamiltonian in a 
frame rotating at the frequency uJi is 



H 



-^ig{a^(T-aa^)^iQ{t){a-a^). (2) 



Here, Ac and are the detuning of the cavity and the 
QD resonance from the laser frequency; is the maxi- 
mum laser strength and p{t) is proportional to the enve- 
lope of the laser electric field. The dynamics of the lossy 
system is determined by using the Master equation 



dp 
'dt 



-i[H, p]^2KC[a]^2jC[a] 



(3) 



where p is the density matrix of the coupled QD-cavity 
system; hi is the cavity field decay rate and 7 is the dipole 
spontaneous emission rate. jC[D] is the Lindblad operator 
corresponding to a collapse operator D. This is used to 
model the incoherent decays and is given by: 

jC[D] = DpD^ - \d^Dp - ^pD'^D (4) 

The Master equation is solved using numerical integra- 
tion routines provided in quantum optics toolbox, trun- 
cating the photon states to 20 photons [8]. This method 
is completely quantum mechanical, and no approxima- 
tion (other than the standard Born-Markov approxima- 
tion and truncation of Fock state basis) is made. 

A semi-classical description of the coupled system [9] 
can be derived by using the relation 



d{D) 
dt 



Tr 



D 



dp 
'dt 



(5) 



valid for any operator D. The mean field dynamical equa- 
tions for the coupled QD-cavity system can then be writ- 
ten as 
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d(a) 
dt 

dt 
d{a.) 
dt 



= -K{a} +g{a) - y^n{t) (6) 

= -7 (<^> + 9 (acTz) (7) 
= -2j{{(T,) + l)-2g{{a^(T) + {a(T^)) (8) 
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FIG. 1: (color online) (a) The schematic of the coupled QD- 
cavity system. It is driven by a laser pulse, and the cavity 
output is monitored, (b) The cavity transmission calculated 
by three different models: the quantum optical (red), semi- 
classical linear (blue) and non-linear (black) model at low 
(Qo/27r = 1 GHz) peak intensity of the driving pulse. All 
three models match quite well. The input pulse is also shown 
(green dashed line). The oscillation in the cavity output is 
due to Rabi oscillation of the photon between the QD and 
the cavity. Inset shows the cavity transmission spectrum in 
presence and in absence of the strongly coupled QD. The split 
resonances are separated approximately by twice the coherent 
dot-cavity interaction strength g. The spectral shape of laser 
pulses with pulse-length 5 ps (blue dashed line) and 40 ps 
(green dashed line) is also shown. Parameters used for the 
simulations are g/27T — 25 GHz, hi/27T — 29 GHz and j/27t = 
1 GHz. 



where, = |e) (e| — \g) {g\. We note that this set of 
equations is not complete and to exactly solve this, we 
need to find the equations describing all the other higher 
order moments, namely (aa^) and (aa^). However, in 
the low excitation regime (no more than 1 photon in the 
system) the QD will remain mostly in its ground state 
and we can approximate (cr^) ^ —1 and replace (aaz) = 
— (a). The resulting set of equations 



d(a) 
dt 

d(a) 
dt 



-f^{a)^g{(j) - ^/f^n{t) (9) 
-j{a)-g{a) (10) 



is identical to the set of equations describing the dy- 
namics of two coupled linear classical oscillators (see 
Appendix A). Although this approximation neglects the 
nonlinear nature of the QD, it matches the actual output 
quantitatively at low excitation power. Unfortunately, 
with increasing drive intensities this model fails com- 
pletely, as the approximation (a^) ^ — 1 becomes invalid. 
For sufficiently high drive intensities though (a^) 
and equation (10) simplifies to 



d(a) 
dt 



-7 (a) 



(11) 



Alternatively, we can retain the dynamics of the a^, 
while making the set of equations (6), (7), and (8) com- 
plete by using the approximations (aa^) ~ (a) (a^) and 



(aV) ^ (a^) (a) [10]. While this approach neglects the 
coherence of the system while analyzing the mean-field 
dynamical equations, the nonlinear behavior of the QD 
is taken into account. 

The temporal cavity outputs at low excitation power 
match very well for the three different models (Fig. lb). 
For the numerical simulation, we used a Gaussian pulse 
with full-width half-maximum (FWHM) of 5 ps as we 
want to drive the dot-cavity system with a pulse having 
bandwidth higher than the coupled system (as shown in 
the inset of Fig. lb). An oscillation in the cavity output 
is observed. This oscillation is due to the coherent Rabi 
oscillation of the photons between the QD and the cavity. 
Note that further oscillations are quenched by the decay 
of the cavity field. 

To intuitively understand the origin of the oscillation, 
we analytically solve the linear semi-classical equations 
to find the eigenvalues of the lossy coupled system as 



±Jg^^-{S-z{^-j)r (12) 



where S = Ac — is the dot-cavity detuning. When 
the real part of the expression under the square root is 
positive, the system is in strong coupling regime and a 
split resonance appears in the cavity transmission spec- 
trum (inset of Fig. 1 b). Without a coupled QD, a 
single Lorentzian peak is observed in the cavity trans- 
mission. The two peaks are entangled states of the QD 
and the cavity, known as polaritons. When the cavity 
is driven with a short pulse, having bandwidth more 
than the coupled system, the cavity output is modulated 
at the frequency difference between the polaritons. i.e., 

2^g^ + l{5-i{K-y)f. 

Although the non-linear semi-classical model allows 
QD saturation, it neglects the quantum mechanical co- 
herence between the QD and the cavity. Fig. 2 com- 
pares the semiclassical and quantum optical simulations 
of the coupled dot-cavity system. We find that the results 
match well both at low (when the QD excited state pop- 
ulation is almost zero and az ^ —1) and high intensities 
of the drive (when the QD is saturated and ~ 0). As 
expected, the nonlinear semi-classical approach will de- 
viate for intermediate intensities. We plot the quantity 
((aV) — {a^){(7))/ftQ integrated over time as a function 
of the driving strength in the inset of Fig. 2. This 
quantity is zero in absence of any coherence. We observe 
that this quantity is smaller for both low and high ex- 
citation power, compared to the value at intermediate 
excitations. Note that the onset of increase in the higher 
excitation power is due to numerical errors caused by the 
truncated Fock state basis. 

Finally, we analyze the dependence of the temporal 
cavity output as a function of four quantities: dot-cavity 
detuning (5, the dot-cavity coupling rate ^, the cavity 
field decay rate and pure QD dephasing rate jd- We 
observe an increase in oscillation frequency (decrease in 
the time interval between the two peaks) when we in- 
crease g (Fig. 3a). This is consistent with the oscillation 
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FIG. 2: (color online) (a) Comparison between the temporal 
cavity transmission obtained via quantum optical (red dashed 
line) and the semi-classical non-linear (black solid line) mod- 
els. The cavity transmission is normalized by the maximum 
cavity transmission and plots are vertically offset for clar- 
ity. The two models match quite well at low and high driving 
power, but at intermediate power, they differ. Inset shows the 
coherence calculated as ((aV) — {a^){a))/Qo integrated over 
time as a function of the driving strength Qq . We observe that 
quantity increases in the intermediate driving power. Param- 
eters used for the simulations: g/27T = 25 GHz and Hi/27T = 29 
GHz. 



period as predicted by the simple linear analysis. At the 
same time, the oscillation period depends only weakly on 

(Fig. 3b). We note that an increasing cavity output 
with increasing cavity decay rate hz. This is due to the 
increasing overlap between the input pulse and the cav- 
ity spectrum. The oscillation frequency increases with 
increasing detuning between the dot and the cavity and 
when the QD is detuned too far from the cavity, the oscil- 
lation almost disappears. This is expected, as with large 
enough detuning the input pulse is not affected by the 
QD (Fig. 3c). An important quantity in solid-state cav- 
ity QED is pure QD dephasing, which destroys the coher- 
ence of the system, without affecting any population of 
the quantum dot states. The effect of pure QD dephasing 
can be incorporated by adding a term 2'yd^{o'^o-) in the 
Master equation [11], where is the pure QD dephasing 
rate. Fig. 3d shows the cavity output as a function of 
pure QD dephasing rate and we observe that the oscil- 
lation eventually disappears with increasing dephasing. 

Finally, we analyze the dependence of the cavity trans- 
mission on the pulse duration (Fig. 4). The pulse dura- 
tion is changed from 5 ps to 50 ps. We observe that 
oscillation frequency is decreasing with increasing pulse 
length. This can be explained by the reduced overlap 
between the input pulse and the coupled dot-cavity sys- 
tem with reduction in pulse bandwidth. In other words. 




20 40 
Time (ps) 



FIG. 3: (color online) The temporal cavity output obtained 
from the full quantum optical simulation as a function of (a) 
the dot-cavity coupling strength g (here k,/27v = 20 GHz; 6 = 
and 7(i = 0), (b) the cavity field decay rate k (here g/27v = 
20 GHz; ^ = and 7d = 0), (c) the dot cavity detuning S (here 
g/27v = k/27v = 20 GHz and 7^ = 0) and (d) the pure QD 
dephasing rate 7cz (here g/27T = hi/27T = 20 GHz and ^ = 0). 
For all the simulations a low excitation power {Qo/27i — 2) is 
assumed. 



a long pulse does not have sufficient bandwidth to ex- 
cite both the polaritons (as shown in the inset of Fig. 
1 b) and the oscillation frequency in the cavity output 
deviates more from 2g. To test the validity of our nu- 
merical simulations, we experimentally probed a strongly 
coupled QD-cavity system. A cross-polarized reflectiv- 
ity setup was used to obtain the transmission of light 
through the coupled system and the cavity transmission 
was monitored with a Hamamatsu streak camera. De- 
tails of the fabrication and the experimental setup can 
be found in Ref. [12] with the experimental parameters 
of the probed dot-cavity system being g/27r = 25 GHz 
and tz/27r = 29 GHz [2]. Unfortunately, we did not ob- 
serve the predicted oscillations in the initial experiments 
measuring the transmission of 5 ps pulses through the 
cavity. This was most likely caused by the limited time 
resolution of our detector. Subsequently the experiment 
was performed with a longer pulse (40 ps FWHM). A 
long pulse does not have sufficient bandwidth to excite 
both the polaritons (as shown in the inset of Fig. 1 b) 
and the oscillation frequency in the cavity transmission 
is different from the value 2^, as shown in Fig. 4. Fig. 5 
a,b,c show the experimentally obtained cavity output for 
three different excitation powers. The experimental data 
match qualitatively the predictions from the numerical 
simulation and clear oscillation is observed in the cavity 
output. This oscillation disappears with increasing laser 
power, as expected from the QD saturation. The oscil- 
lation period is estimated to be 25 GHz, corresponding 
to time difference of 39 ps between the two peaks. We 
note that the numerically obtained plots in Fig. 4 are 
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FIG. 4: (color online) The normalized cavity transmission for 
different pulse duration. The pulse duration is changed from 
5 ps to 50 ps. We observe oscillation in the cavity output, 
although the oscillation frequency decreases with increasing 
pulse-width. This can be explained by the reduced overlap 
between the pulse and the coupled dot-cavity system in fre- 
quency domain. 



done with very small excitation power. However, the ex- 
periment cannot be performed with such low excitation 
power as the detected signal is too low. Hence, in the ex- 
periment, the coupled system is driven close to the QD 
saturation, and the oscillations are much less visible. 

In summary, we have analyzed the nonlinear temporal 
dynamics of a strongly coupled QD-cavity system driven 
by a short laser pulse. We showed that this quantum 
optical system behaves similar to two coupled classical 
linear oscillators when the system is driven with a weak 
pulse and that a signature of the vacuum Rabi oscillations 
can be observed in the time resolved cavity transmission. 
For high intensity pulse these oscillations die down due to 
saturation of the QD. We provided a semi-classical non- 
linear model and showed that in the actual dynamics 
the role of quantum coherence is important. Lastly, we 
presented experimental evidence of those oscillations in 
the cavity output. 
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Appendix A: Dynamics of two classical coupled 
linear oscillators 

The dynamics of two classical coupled oscillators, with 
resonance frequency ujo and decay rates Fi and are 
governed by 



Pi 



dxi 



- uJqXi + G{xi — X2) = ^{t)e 



iojot 



and 



X2 



dt^ 



dX2 

'~dt 



-UjIx2 



G{x2 - xi) 0, 



(Al) 



(A2) 



where G denotes the coupling strength between the os- 
cillators. One of the oscillators is driven resonantly with 
driving strength l^(t), as the cavity is driven by a laser. 
We assume the solution of the form xi{t) = Xi(t)e*^°^ 
and X2{t) = X2(t)e*'^o^ where Xi{t) and X2{t) are slowly 
varying envelopes of the actual oscillator outputs. Then 
we can write 



dxi 
~dt 
d^xi 
~dt^ 



' dX 

iujoXie'"^^' + ( ^e^^o' 
dt 



dX 

dt ° 



For X2 we can find similar equations. Using slowly 
varying envelope approximation << iuj^Xi and 

^^^^ << zujo^^^ujqXi)^ we remove the bracketed terms 
and get the following equations for the un-driven coupled 
oscillator system: 



dXi 
~dt~ 



£1 
2 



G 
2iujQ 



G 



X2 + n{t) (A3) 



and 



dX2 
dt 



r2 

2 



G 

2iu)r 



X2 



G 
2iu)o 



Xi 



(A4) 
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FIG. 5: (color online) Experimentally measured time resolved transmission of 40ps pulses through a strongly coupled dot-cavity 
system for three different powers (averaged over the pulse repetition period): (a)O.l nW; (b) 0.23 nW; (c) 1 nW. The powers 
are measured in front of the objective lens in the confocal microscopy setup. For this specific system cavity field decay rate 
kI2i^ — 29 GHz and coherent dot-cavity coupling strength gjlix — 25 GHz. Clear oscillations are observed in the cavity 
transmission, consistent with the the theoretical predictions. We also observe decreasing oscillation with increasing laser power 
due to QD saturation. 
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